Photon counting with loop detector 



Konrad Banaszek and Ian A. Walmsley 

Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom 

(Dated: February 1, 2008) 

We propose a design for a photon counting detector capable of resolving multiphoton events. The basic 
element of the setup is a fiber loop, which traps the input field with the help of a fast electrooptic switch. A 
single weakly coupled avalanche photodiode is used to detect small portions of the signal field extracted from the 
loop. We analyze the response of the loop detector to an arbitrary input field, and discuss both the reconstruction 
of the photon number distribution of an unknown field from the count statistics measured in the setup, and the 
application of the detector in conditional state preparation. 
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PACS numbers: 040.5570, 270.5290 

Photon counting is an important method of detecting light 
with applications in many fields including spectroscopy, at- 
mospheric physics, and quantum information processing. De- 
tection of single photons necessarily involves a gain mecha- 
nism in order to generate a macroscopic number of photoelec- 
trons from the absorption of a single energy quantum from the 
incident electromagnetic radiation. Among the most popular 
of detectors with this characteristic are avalanche photodiodes 
operated in the Geiger mode. Compared to standard pho- 
tomultipliers, they have higher quantum efficiency, are more 
stable and are more robust with respect to external conditions. 

Geiger-mode avalanche photodiodes have one drawback: 
the breakdown current is almost completely independent of 
the number of absorbed photons. Consequently, it is not pos- 
sible to determine the number of photons incident on the de- 
tector on a time scale short compared to the detector response 
time. This has serious implications for several applications. 
For example, it critically affects the performance of practical 
quantum cryptography systems. [||| It also impairs the fidelity 
with which entangled states may be prepared, [Q| as well as 
the conditional detection required for linear optical quantum 
computing. [Q| 

In this paper, we propose a photon counting setup which 
is capable of resolving multiphoton events. The basic idea is 
to split the input signal into separate small pieces which are 
expected to contain less than one photon, and therefore can 
be detected with an avalanche photodiode without losing in- 
formation on the photon number. Our method is based on 
current fiber optic technology and has the important advan- 
tage that it requires only a single avalanche photodiode. This 
contrasts with previous proposals that involve the splitting the 
signal input, which require a large array of photodiodes. [Qj 
We note that multiphoton absorption events can be in princi- 
ple resolved by more exotic photon-counting detectors, such 
as solid state photomultipliers.[^] These, however, still remain 
at the developmental stage, and must be operated at liquid he- 
lium temperatures. 

The proposed loop detector is depicted schematically in 
Fig. [l|. The basic element of the setup is a fiber loop used 
as a storage cavity. The light pulse whose photon statistics 
are to be measured enters the loop through a fast electrooptic 
switch S. Immediately after injection of the pulse the state of 
the switch is flipped, which closes the loop and thus traps the 
input pulse. In each circulation of the loop, a small portion of 



the pulse is extracted by a weak static coupler C and incident 
on avalanche photodiode (APD) operated in the Geiger mode. 
The duration of a single roundtrip in the loop is longer than the 
dead time of the APD, so that after each roundrip the APD is 
ready for a new detection event. The quantity recorded in the 
measurement is the total number of counts on the photodiode. 

The main result of this paper is to show that the statistics 
of counts on the APD are related to the photon distribution of 
the input pulse. We will first calculate the response of the de- 
tector to a coherent field. This result can be generalized to an 
arbitrary input field by averaging with an appropriate proba- 
bility distribution for the field amplitude. In the most general 
case this probability distribution is given by Glauber's P rep- 
resentation for the input field. Next, we will use this result 
to derive a method for reconstructing the photon number dis- 
tribution from the count statistics for a completely unknown 
input field. 

We begin by calculating the probability distribution pi(k) 
of obtaining k counts on the photodiode assuming that the 
coherent pulse after injection into the loop had intensity T. 
In the derivation, it will be convenient to use the generating 
function: 



k=0 



(1) 



There are several parameters of the setup relevant to our cal- 
culation. Let t r be the power transmission of the loop for a 
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FIG. 1: The proposed construction of the loop detector. S, electroop- 
tic switch; C, weak coupler; APD, avalanche photodiode. 
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single pulse roundtrip, including losses at all the optical ele- 
ments. Further, let t c be the fraction of the light intensity ex- 
tracted by the coupler C to the avalanche photodiode. We as- 
sume that the photodiode is characterized by the quantum ef- 
ficiency r), and that the probability of registering a dark count 
is pd- We will discuss realistic values for all these parame- 
ters later, though our analysis is valid for an arbitrary set of 
parameters satisfying certain approximations. 

The assumption of coherent input substantially simplifies 
the calculations, since the portion of the pulse directed to the 
APD in each roundtrip is uncorrected with the field remain- 
ing in the loop. Consequently, the counts registered by the 
avalanche photodiode in each roundtrip are statistically inde- 
pendent. In the zth circulation of the loop, the intensity of the 
pulse is Of this the light intensity extracted to the pho- 

todiode is t c t % ~ l I. Taking the efficiency of the avalanche pho- 
todiode to be r], the probability of a count in the ith roundtrip 
is 1 — (1 — pd) exp(—rjt c t' l r ~ 1 T), including the possibility of a 
dark count with probability pd- For an input pulse in a coher- 
ent state the counts on the APD are statistically uncorrelated 
on each circulation, so that the generating function p%{z) is 
simply given by a product of generating functions describing 
detection events for each round trip: 

L 

Vt{z) = lit 2 + C 1 - ~ eM-VtctlT 1 !)}. (2) 

i=l 

Here L is the total number of circulations of the loop. 

The expression derived in Eq. (||) describes the most gen- 
eral situation of arbitrary input intensity. An important lim- 
iting case is when the probability of extracting two or more 
photons from the loop in a single roundtrip via the coupler 
C is negligible. Then the logarithm of the generating func- 
tion lnpx(z) can be expanded up to terms linear in t c I or pd, 
assuming that both these parameters are much smaller than 
one. Further, when the number of roundtrips is large enough 
to allow all the input signal photons to leak from the loop, we 
obtain from the logarithmic expansion the following approxi- 
mate expression for the generating function: 



pi(z) w exp 



(*-l) 
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(3) 



This formula describes simply the standard Poissonian statis- 
tics normally associated with coherent radiation, with the av- 
erage number of counts equal to rjt c X/(l — t r ) + Lpd- The 
first term here describes the average number of counts gen- 
erated by the input light, and the second term corresponds to 
dark counts. In this regime, the response of the loop detector 
is described by a single effective efficiency parameter: 



'/eff = I] 



1 - t r 



(4) 



In the ideal limit all the optical elements are lossless, and the 
only source of attenuation for the pulse trapped inside the loop 
is extraction of the photons to the avalanche photodiode. Then 
t c = 1 — t r , and the effective efficiency of the loop detector 
approaches that of the avalanche photodiode, but with the ad- 
vantage of having photon number resolution. In a realistic 



case, there are always excess losses in the coupler or the elec- 
troopic switch, which decreases the effective efficiency of the 
loop detector below the level of the APD itself. 

We now turn to the problem of reconstructing the photon 
number distribution of an unknown input field from the statis- 
tics of counts measured using the loop detector. Let g(n) be 
the photon number distribution of the pulse immediately after 
injection into the fiber loop. The probability p(k) of observing 
k photocounts is given by the linear combination: 



Pi k ) ^^2w(k\n)g(n) 



(5) 



where w(k\n) are conditional probabilities of registering k 
counts on the photodiode given exactly n photons injected 
into the loop. These conditional probabilities are defined ex- 
clusively by the parameters of the setup. The explicit form of 
these conditional probabilities can be found from our previ- 
ous calculation of the generating function for a coherent input. 
For a coherent input pulse, the photon number distribution is 
Poissonian, g(n) = e~ x T n /nl. Inserting this into Eq. (^), 
multiplying both its sides by z k and performing the summa- 
tion over k yields: 



oo oo 



e I Pi{z) = ^2^2 w(k\n) 



z k 2" 



(6) 
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where px(z) is defined in Eq. (||). Thus the conditional proba- 
bilities w(k\n) can be identified from the generating function 
e I px(z) for a coherent input by expanding it into a double 
power series in the parameters T and z. This procedure can 
be easily performed using any of the standard computer pro- 
grams for symbolic algebraic calculations. Given the explicit 
form of the conditional probabilities w(k\n) and the measured 
count statistics p(k), one can then use one of the standard 
methods for solving linear systems to reconstruct the photon 
number distribution g(n). In Fig. |we show the results of a 
Monte Carlo simulation of the loop detector operation, using 
the singular value decomposition method to reconstruct the 
photon number distribution. The detector is able to provide 
very good estimates of the photon number distribution of both 
classical and nonclassical inputs. The accuracy of the recon- 
struction can be estimated using standard statistical tools. [Q] 

In a number of quantum information processing 
protocols, ^ [i]] determination of the photon number has 
to be made on the single-shot basis in order to infer non- 
destructively the state of another quantum system that is 
entangled with the detected radiation. In such an application, 
the performance of the detector can be characterized by the 
confidence: [|[ 

w(k\k)g(k) 



(7) 



J2 n w(k\n)g(n) 

describing a posteriori probability that a fc-click event has 
been triggered by the fc -photon component of the input field 
state. Here g{n) is the a priori photon number distribution 
of the input field, and w(k\n) are the conditional probabili- 
ties calculated before. In practice, the confidence can be opti- 
mized by tuning the splitting ratio of the coupler C used in the 
construction of the loop detector. 
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FIG. 2: Monte Carlo simulation of photon counting using a loop de- 
tector for a coherent state with the average photon number equal to 
one (left), and a single-photon Fock state (right). The top graphs 
depict simulated count statistics based on N — 10 s repetitions of 
the mesurement, assuming t r = 72%, t c = 20%, r\ = 80%, 
Pd = 0. The dots in bottom graphs present the photon number distri- 
butions reconstructed with the help of the singular- value decomposi- 
tion method, compared with the exact distributions shown with bars. 



An important question is whether realistic values for the pa- 
rameters of the loop detector enable it to function in the man- 
ner shown in Fig. g. In the visible region, silicon avalanche 
photodiodes can exhibit efficiencies greater 75% for wave- 
lengths near 700 nm.[|[ ||] The dark count rates for these de- 
vices can be as low as 25 Hz, which is negligible. This rate 
may be further reduced by operating the photodiode in the 
gated mode. The APD parameters become less favorable be- 
yound the 1 wavelength region, where germanium or In- 
GaAs materials are used. [|To[, pTTj , |i"2"| ] As with single photon 
detection, the quantum efficiency of the photodiode imposes 



the upper limit on the performance of the loop detector. A 
second important factor is the excess losses inside the fiber 
loop, which attenuate the pulse in each roundtrip in addition 
to the fraction extracted by the coupler C. These excess losses 
lead to two opposing constraints on the extraction of the light 
from the loop to the photodiode, parameterized by the cou- 
pling losses t c . On one hand, we have seen that in order to 
recover the complete photon number distribution, t c should be 
as small as possible. On the other hand, if too small a power is 
extracted from the loop, then a large number of roundtrips are 
required to measure the count statistics. Consequently, sub- 
stantial excess losses in the loop would mean that most of the 
photons would disappear inside the loop before ever reaching 
the photodiode, and would give no signal at all. The optimal 
value of t c is some intermediate value which allows one to de- 
tect a substantial fraction of photons in the input signal and at 
the same time provides a sufficiently large number of multi- 
photon events. A good test of whether both these conditions 
are satisfied is to check the singularity of the matrix composed 
of the conditional probabilities p(k\n) for the range of photon 
numbers which are expected in the input signal. If the ma- 
trix is not singular, then the count statistics provides sufficient 
information to reconstruct the complete input photon number 
distribution from the experimental data. 

In conclusion, we have proposed and analyzed a design for 
a photon counting detector capable of resolving multiphoton 
detection events using standard laboratory technology. The 
setup uses commonplace fiber-optic components and does not 
require extreme operating conditions. 
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